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An atomistic insight into potential barrier formation and band bending at the interface between 
a metal and an n-type semiconductor is achieved by ab initio simulations and model analysis of a 
prototype Schottky diode, i.e., niobium doped rutile titania in contact with gold (Au/Nb:Ti02). 

The local Schottky barrier height is found to vary between 0 and 1.26 eV depending on the position 
of the dopant. The band bending is caused by a dopant induced dipole field between the interface 
and the dopant site, whereas the pristine Au/Ti02 interface does not show any band bending. 

These findings open the possibility for atomic scale optimization of the Schottky barrier and light 
harvesting in metal-semiconductor nanostructures. 


The presence of a Schottky barrier (SB) between a 
semiconductor and a metal is of paramount importance 
to numerous application fields, including electronics [ljj, 
photovoltaics |2|, y], photocatalysis [J-lg and gas sen¬ 
sors 0-0. Schottky barrier physics has been a subject of 
intense investigation for several decades, but has recently 
received renewed substantial attention in two areas: i) 
the emergence of novel Schottky devices in plasmonics 
for photocurrent generation, photo detection and solar 
light harvesting aim Ei; and ii) the development of 
quantum-scale metal-semiconductor structures, pushed 
by the ever present need to further minimize and opti¬ 
mize electronic devices fl2l - [l7j . Continued development 
of these areas could be greatly facilitated by an atomistic 
understanding of SB-based processes. 


The quantum transmission of electrons or holes across 
the SB is determined by two quantities: the barrier height 
and, more importantly, the decay length of the band 
bending. Together, these quantities determine the prob¬ 
ability of transmission and the energy distribution of hot 
carriers across the metal-semiconductor interface. The 
conventional SB model assumes a uniform charge deple¬ 
tion region on the semiconductor side and a charge accu¬ 
mulation layer localized at the interface 00 , result¬ 
ing in a parabolic bending of the semiconductor bands. 
The decay length of the band bending has been believed 
to be on the order of 10 nanometers for typical dopant 
concentration (10 19 cm" 3 ). However, recent advances in 
nanotechnology 00 Elj] . which has made it possi¬ 
ble to control and characterize the SB at the nanome¬ 
ter scale, have revealed important deviations from pre¬ 
dictions made from the homogeneous Schottky barrier 
height (SBH) model 14|. Instead, the results, which de¬ 
pend on materials properties, dopant compositions and 
concentrations, have been qualitatively interpreted in the 
inhomogeneous SBH model 0, 0, although the cur¬ 
rent lack of a complete atomistic picture severely limits 
the possibility to achieve a quantitative understanding 
of SB formation. In particular, conventional models all 
assume a uniform dopant distribution, and it remains 


unclear how the composition and atomic structure of the 
semiconductor affect the electronic structure, i.e., barrier 
height and band bending, on the atomic level. 

Here, we report first-principles calculations of a pro¬ 
totype Schottky diode, i.e., the Au/Ti02 interface. A 
substitutional niobium dopant (Nbri) was introduced 
to model discrete defects in the n-type semiconductor. 
The concentration and spatial distribution of the NbTi 
can be controlled experimentally 24H27)] and Nb dopants 
are known to induce small lattice relaxations [28fl. which 
makes this particular system suitable for atomistic sim¬ 
ulations. We found that the pristine Au/Ti02 interface 
has a relatively large barrier height, but shows no band 
bending. The band bending and decay length is instead 
determined by the precise locations of the dopant. Our 
results also show that the band bending is inhomogeneous 
and highly localized to the defect region. We calculated 
the dopant position dependent barrier height and show 
that it can be qualitatively understood by the deep level 
(DL) barrier model 0|. Our results reveal the origin 
and nature of inhomogeneity of the SBH and shed light 
on the mechanisms of electron transmission across the 
metal-semiconductor interfaces. 


All the calculation in this work was done using the 
Vienna Ab initio Simulation Packages (VASP) 0-33 
with the projector augmented wave (PAW) method [34, 
|35| and the PBE 0,0] exchange-correlation functional 
in the generalized gradient approximation (GGA). The 
Coulomb correlation of the Ti 3d orbitals were treated in 
the GGA+U scheme with an effective on-site Coulomb 
repulsion U=10 eV, which was calibrated to the band 
gap and energetic position of the defect states. This U 
parameter results in a 3.28 eV indirect band gap at M—T 
and a 3.36 eV direct band gap at T for bulk rutile Ti02- 
The Nb defect states are located below the conduction 
band minimum (CBM), (see Supplemental Material 0] 
for calculation details and the choice of U parameters). 

In order to explore the effect of chemical composition 
on the SB, a large number of different atomic struc¬ 
tures were calculated based on the lattice alignment 
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and orientation of Au nanoparticles on thin film rutile 
(110) TiC >2 [39|]. In epitaxial growth, the Au[110] close- 
packed direction is always found to be parallel to the 
TiO 2 [ 001 ] direction with a lattice mismatch as small as 
0.4%. Different epitaxial layers are obtained by rotating 
around the [110] axis. In the following, we focus on the 
Au(112)/TiO2(H0) interface, which has been observed 
after high temperature deposition 0- In this case, one 
Au atom is located on top of Ti 5c in the rutile TiO2(110) 
surface (Fig. QJa)) [40]. We modeled the interface in a 
slab geometry using a 12 O-Ti-O trilayers (40 A thick) 
in contact with 3 layers (5 A thick) of Au(112). The 
periodic slabs were separated by a 15 A thick vacuum 
region (Fig. Oja)), and dipole correction was used to al¬ 
leviate the image interactions. The work function of the 
Au slab (^m) and the electron affinity of the TiC >2 slab 
(x«) were calculated to be 5.13 eV and 4.27 eV, which is 
in good agreement with experimental values (5.2 eV 41 ] 
and 4.3 eV >42], respectively). The Nb-dopant was in¬ 
troduced by replacing one of the Ti atoms in the lattice 
(Nbri)- In addition, the effect of changing the concen¬ 
tration of the Nb-dopants was checked by increasing the 
surface unitcell to (1 x 2) and (2 x 3) using a 5-layer TiC >2 
slab. 

The electronic structure of the pristine interface is 
shown in Fig. [T| The lower panels show the band struc¬ 
tures of the isolated Au and TiC >2 slabs (left) and the in¬ 
terface formed upon contact (right). The interface band 
structure reveals a rigid upward shift of about 0.5 eV for 
both the valence and conduction bands of TiC> 2 , while 
the Fermi level and Au bands change by less than 0.1 eV. 
This difference is caused by charge polarization, as shown 
in panel (b). We analyzed this effect using the contact 
induced charge density p and doping induced charge 
density pf nd defined as: 


Pind = PMS - (pm + Ps), 

Pind = PD-MS ~ (PMS ~ PTi + PNb), 


(1) 

( 2 ) 


tCd = 


P™d /D ( z - z o)dz. 


(3) 


Here, pMS and pd-ms are the charge densities of the un¬ 
doped and doped Au/TiC >2 interfaces, respectively, pM 
and ps are those of the isolated metal and semiconduc¬ 
tor slabs, while pxi and pNb are the charge densities of 
the free atoms. The areal density of the interface induced 
dipole p \™ d was calculated by integration over an inter¬ 
face range of 5 A and found to be = — 0.014.DA 
The negative sign means that the dipole points from Ti02 
to Au (Fig. [ljb)). The charge transfer between Au and 
undoped TiC >2 is found to be less than 0.05 e/cell, as 
shown by a Bader analysis. The induced dipole results 
in a SBH of 1.26 eV, which is 0.4 eV larger than the 
difference between the work function and the electron 
affinity of the isolated Au and TiC >2 slabs. The band gap 
of TiC >2 is 3.3 eV, that is the same as in the bulk calcu¬ 



FIG. 1. The geometry and electronic structure of the pris¬ 
tine interface show that the TiC >2 band shifts up owing to 
the influence of the induced interface dipole, (a) Schematic 
of the Au (112) and Ti02 rutile (110) interface used in the 
calculations. The boxes show the cell employed in the peri¬ 
odic boundary calculations, (b) The plane (parallel to the 
slab) averaged induced charge density. Red and blue indi¬ 
cate charge accumulation and depletion. The induced charge 
density gives rise to a dipole at the interface, (c) The band 
structure of isolated (left) and combined (right) Au and rutile 
TKH slab. The red dashed line is the Fermi level. The side 
columns show the density of states (DOS) of the Au and WO 2 
components. The Ti02 bands shift up by 0.4 eV owing to the 
contact with Au. 


lation. The interface position (z = 0) was chosen at the 
electronic potential maximum between the TiC >2 and Au 
layers. Importantly, the conduction band edge is found 
to be flat throughout the semiconductor slab (Fig. S3). 

The results above show that the pristine Ti02/Au in¬ 
terface does not exhibit any band bending in the semi¬ 
conductor region, and we therefore hypothesized that the 
atomic scale SBH inhomogeneity and band bending in¬ 
stead are caused by the dopant. Figure [5] shows a typical 
case with a Nb-dopant placed two layers beneath the in¬ 
terface. Here the SBH is defined as the energy dif¬ 

ference between the Fermi level ( Ef ) and the conduction 
band minimum (CBM) of the TiC >2 layer at the interface. 
We found that both the potential profile and the SBH 
critically depend on the position of the dopant. When the 
Nbn is located at the contacting layer, the defect state 
density is high enough to pin the Fermi level, which is 
very close to the bottom of the conduction band, as found 
in previous studies j4f|. In this case, the SBH is signifi¬ 
cantly reduced, and the band bending is the same as for 
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FIG. 2. Illustration of Schottky barrier formed at a 
Au/Nb:Ti02 interface. The color map shows the local density 
of states (LDOS), as obtained from the DFT+U calculations, 
along the direction normal to the interface. The Schottky bar¬ 
rier height for n-type doping is defined as the energy 

difference between the Fermi energy (red line) and the con¬ 
duction band minimum (upper purple line) at the interface. 
The Nb-dopant is marked as Nbn in the inserted structure 
plot. The dopant state below the bottom of the conduction 
band traps most of the excess electron donated by Nbn■ The 
band edges of the pristine Ti 02 are also shown (gray dashed 
lines). The band bending is caused by the dopant induced 
charge polarization. 


the pristine case, i.e. it is essentially flat. However, as 
the dopant is moved away from the interface, the dopant 
induced charge is distributed nonuniformly between the 
interface and the dopant, leading to a lowering of the 
potential towards the dopant, i.e., a clear band bending. 
This trend starts from the second layer (dopant-interface 
distance do = 5.9 A) and becomes more prominent as 
the dopant moves to the other end of the slab at 12 th 
layer (do = 39.2 A). Thus, our results clearly show that 
the local dopant confine the SBH and dictate the band 
bending. 

The variation in charge and potential with dopant- 
interface distance are analyzed in Fig. [3] With the Nb- 
dopant positioned between the second layer and the fifth 
layer, the band bending is nearly parabolic. This is 
consistent with the uniform dopant SB model. How¬ 
ever, it becomes essentially linear if the Nb-dopant is 
located beyond the fifth layer, i.e. da > 16 A. The 
magnitude of the induced interface dipole is found to be 
inversely proportional to the dopant-interface distance, 
= 0.066 DA -1 /(dr, - 1.08 A) (Fig. S4), and the 


FIG. 3. The dopant-interface distance (do) dependence of the 
induced charge and band bending, (a) Nbn induced charge 
density. The curves are the plane averaged value and the color 
inserts are the isosurface with isovalue = 10~ 3 A 3 . Yellow 
(+), Green (-). The induced interface dipole is inversely pro¬ 
portional to do. (b) Band bending in Ti02- The points are 
the local conduction band minimum (CBM) in each layer in 
relation to the Fermi energy. The vertical bar indicates the 
position of Nb-dopant in each calculation. The Ti02 band 
shifts up while approaching the interface. 


dipole points from Au to Nb:TiC> 2 , while the induced 
dipole around the dopant points in the opposite direction. 
The charge transfer between the metal and the semicon¬ 
ductor slabs was found to be negligible for the pristine 
interface as well as for the doped cases (0.07 e/cell in case 
of interfacial Nbxi, and less than 0.02 e/cell at other 
doping positions). The excess electron contributed by 
the Nb-dopant is instead mainly distributed on the Nb 
atom and the neighboring Ti atoms along the [001] di¬ 
rection. This localized charge distribution is consistent 
with previous DFT+U calculations 28] and STM exper¬ 


iments [44]. Further Fig. [31(b) shows that the CBM of 
the layer containing the dopant was pinned to the energy 
Ef + C eV, where £ is in the range of 0.2 ~ 0.5 eV and 
varies with do- The same shift was observed for the Ti 
4s semi-core states (Fig. S5). Thus, we can conclude 
that dopant induced charge polarization give rise to a 
dipole field that is mainly responsible for the band bend¬ 
ing. We can further conclude that the SBH and band 
bending are highly inhomogeneous (locally determined) 
and strongly dependent on the dopant-interface distance 
within a range of a few nanometers. 

The atomic picture described above is in sharp con¬ 
trast to the conventional SB picture, which does not in- 
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elude any dopant position dependence as it assumes that 
the band bending is caused by the electrostatic poten¬ 
tial generated by uniform ionized defects in the deple¬ 
tion layer E m The conventional potential profile 
consists of aparabolic term plus an image charge correc¬ 
tion term 1 1811451 1 - The image force (IF) decreases the 
SBH by 


A &I n = 


e 6 N, 


D 


8t r 2 ege| 


- c - k B T) 


(4) 


(a) 1.5 


0.5 



d D (A) 


where e is the charge of the electron, Nd is the dopant 
concentration, £o and Es are the vacuum permittivity 
and relative dielectric constant of semiconductor, respec¬ 
tively, n is the SBH in the absence of image charge 
correction, £ is the energy difference between the CBM 
and the Fermi level, is Boltzmann’s constant and T is 
the temperature. Using the geometry from our first prin¬ 
ciples calculations, i.e. one Nbri in the long (1 x 1 x 12) 
cell, we have Nd = 1.3 x 10 21 cm~ 3 . The depletion layer 
width is found to be 10 A (25 A) with es = 10 (es = 60). 
The SBH reduction is then 0.58 eV for es = 10 (0.15 eV 
for es = 60) (Fig. S6), which deviates from the calcu¬ 
lated SBH reduction (Fig. dj), and does not explain the 
dependence on the dopant position (Fig. [3]). This qual¬ 
itative and quantitative discrepancy between the SBH 
prediction obtained from the uniform dopant Schottky 
model and our results clearly indicate the importance of 
an atomistic description of the interface. 

We now compare the DFT+U results with the alterna¬ 
tive deep level (DL) model [ 29 }, which was developed to 
specifically incorporate the local barrier profile near the 
interface. The DL model assumes point charge donors 
with energy ( Edl ) below the CBM. The SBH reduction 
depends on the charge donor-interface distance (dpe) 
and the areal dopant density ( ctdl) according to [29( 


A$ 


DL _ 
B ,n 


e 2 aa DL (l - f DL ) 
£0 £S 


(5) 


where 

Idl = 1 1 + exp 



—e 2 d D L<TDL (1-/di)/(eoEs)-®M 

FeTT 



is the Fermi-Dirac distribution function of the DL states. 
We extracted the relevant parameters from the DFT+U 
calculations (Fig. [31): n = 1.26 eV is the SBH without 

DL states, a = 5 A is the characteristic width of the in¬ 
terface, Edl = 0.2 eV is the DL energy below the CBM, 
and we set es = 10 and T = 1000 K. The resulting SBH 
obtained from the DL model is plotted in Fig. [4] and is 
found to be in good agreement with our first principles 
calculations, except for the regions closest to the bound¬ 
aries. To pin down the origin of this discrepancy, we 
also performed calculations for a case where the atomic 
structure of the Nb-doped slab were fixed at the pristine 
interface structure position (DFT+U (fix geo.) in Fig. [4}. 
The DL model then captures the main dopant dependent 


FIG. 4. Schottky barrier height reduction (A &B,n) as a func¬ 
tion of (a) dopant-interface distance (do) and (b) dopant areal 
concentration (od). The histogram show the DFT+U results 
obtained with full geometry optimization. The circles are 
calculations with atomic coordinates fixed at the pristine in¬ 
terface geometry. The solid lines are calculated using the DL 
model (Eq. [5j while the dashed lines are obtained from the 
uniform dopant Schottky model including an image force (IF) 
correction (Eq. [4j. 


SBH feature extremely well, the only exception is the case 
when the dopant is located at the interfacial layer where 
the DL model is not applicable. 

In conclusion, first principles calculations of Au/TiCU 
interfaces show that the SBH is tuned by interface dipoles 
induced by contact and dopant. The local barrier pro¬ 
file, i.e. the band bending, shift almost linearly between 
the interface and the dopant layer. The barrier width 
is a-few-layer thick and depends on the dopant position. 
The reported experimental SBH for Au/TiCU interfaces 
lies in the range 0.9~1.2 eV [§, _16, 147) . Given that these 
measurements are macroscopically averaged, our calcu¬ 
lations are in very good agreement with experiment. In 
contrast to the uniform dopant Schottky model, the DL 
model is able to account for the SBH reduction varia¬ 
tion with dopant position. The conclusion and overall 
picture emerging from the present study should be gen¬ 
erally applicable and highly relevant also to other metal- 
semiconductor systems. As such, they can serve as a 
basis and reference to further studies of internal electron 
emission and hot-carrier transport across metal semicon¬ 
ductor interfaces. In light of the rapid development of 
layer-controlled molecular beam epitaxy (48} , the results 
open up the possibility for atomic scale engineering and 
optimization of novel SB-based devices. 
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COMPUTATION METHOD 

In the density functional theory (DFT) calculations, 
the plane wave cut-off energy was set to 400 eV. A 
Monkhorst-Pack k-point mesh of 5 x 11 x 1 centered at the 
r point was used for the Au/TiCb rutile (110) interface 
supercell. A more condensed k-sampling (9 x 21 x 1) was 
adopted when extracting the density of states (DOS). 
The DFT+U calculations were performed using Du- 
darevs scheme [l|. The geometry optimizations were 
done on the GGA level. All of the atoms were relaxed 
in the geometry optimization processes, and were con¬ 
sidered to be relaxed when the maximum forces on each 
atom was smaller than 0.01 eV/A. The calculated lattice 
constants of Ti02 (rutile) and Au (fee) are a = 4.64 A, 
c = 2.96 A, u = 0.305 and a = 4.17 A, respectively, in 
good agreement with experimental values. 

U DETERMINATION 

In the DFT+U calculations, the on-site Coulomb re¬ 
pulsion U parameter was carefully calibrated by calculat¬ 
ing the bandgap of pristine rutile Ti 02 and comparing 
with results from both Go Wo calculations and experi¬ 
mental measurements. The calculations were performed 
using a cell with 6 atoms and a 9 x 9 x 11 Monkhorst- 
Pack k-point sampling. The parameter was determined 
to be U(Ti3d) = 10 eV (Table ISTI) . Furthermore, with 
this value the electron affinity of 12-layer Ti 02 rutile slab 
was found to be 4.27 eV, whichjdso agrees well with the 
experimental value of 4.33 eV fH- 

TABLE SI. The band gap ( E gap ) in bulk rutile TiCU and 
electron affinity (y s ) of the 12-layer rutile (110) slab, (units 
in eV) 


redirect rpindirect ~ 

-‘-'gap -‘-'gap Xs 

(r ->■ r) (m r) 



0 

1.88 

- 



2 

2.12 

- 


DFT+U 

4 

2.39 

- 


U{Ti3d) = 

6 

2.69 

2.67 



8 

2.99 

2.95 



10 

3.31 

3.25 

4.27 

Go Wo 


3.41 

3.34 

4.33 0 

Exp. 
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The U(NbAd) was then determined by keeping the 
U{Ti3d) fixed at 10 eV and calculating the dopant state 
with U(NbAd) varied from 0 eV to 5 eV (Table IS2l) . A 
supercell of 2x3 x 2 along [110], [001] and [110] directions 
was considered. The Nb dopant occupies one of the host 
Ti sites within the supercell Q, thus forms Tii_ x Nba;02 
with x = 0.042, which is equivalent to our 12-layer Ti02 
calculations (a dopant concentration of 1.3 x 10 21 cm -3 
and 4.7 wt%). Higher concentration of Nb-dopant was 
also calculated using the same supercell but with two 
Ti atoms replaced by Nb. With U(NbAd) < 3, the 
charge carrier distribution (Table IS2l) and states at the 
bottom of conduction band (Fig. ISTI) were not sensitive 
to the varied value of the U(NbAd) parameter. When 
U(NbAd) increases from 0 to 3 eV, a dopant state lo¬ 
cates at 1.3 ~ 1.5 eV below the conduction band edge 
and about 0.5 eV below the Fermi level, while the ex¬ 
tra charge was trapped on neighbor Ti atoms along [001] 
direction. The DOS below the conduction band edge 
(CBE) was found compatible with x-ray photoemission 
spectra (XPS) 0 (Fig. 1521) . and the charge transfer to 
the neighboring Ti atoms along [001] direction created 
a polaron in agreement with previous studies jBJ. For 
U(Nb4d)> 4 eV, the dopant state is upshift to just be¬ 
low the conduction band edge, which results in remark¬ 
able changes in the charge distribution (Tabic IS2l) . Based 
on the consideration of the defect state energy distribu¬ 
tion and charge carrier localization, U(NbAd) = 0 was 
sufficient in this work. 


TABLE S2. The space distribution of the extra electron in¬ 
troduced by Nbn on Nb atom and the two nearest neighbor 
Ti atoms along [001] directions. x2 indicates the charge dis¬ 
tributions 011 the two nearest neighbor Ti atoms are equal. 


U(Nb4d) 

UNb 

UTi 

0 

0.101 

0.363 x 2 

1 

0.091 

0.368 x 2 

2 

0.082 

0.374 x 2 

3 

0.073 

0.378 x 2 

4 

0.451 

0.154 x 2 

5 

0.622 

0.064 x 2 
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FIG. SI. The projected density of states (PDOS) of 
Nb x Tii- x C >2 on host ( Ti\- x C> 2 ) and dopant ( Nb ). The 
U(Ti3d) parameter was fixed at 10 eV and the U(Nb4d) var¬ 
ied from 0 to 5 eV. The energy was aligned according to Ti 
4s semi core orbital energy. The vertical red bars mark the 
Fermi energy in each calculations. 



FIG. S2. The calculated density of states (DOS) using 
U(Ti3d) — 10 eV and U(Nb4d) = 0 eV compared with XPS 
experiments Q. The energy was aligned to the conduction 
band minimum. 


STEP-LIKE BARRIER AT PRISTINE 
INTERFACE 


Pristine Au/TiC >2 interface forms step-like barrier as 
shown in Fig. IS3I The electron polarization at the inter¬ 
face forms a plane of induced dipole (Fig. 1(b)). The 
potential in the semiconductor side upshifts by a con¬ 
stant value. Figure [S3] shows the electrostatic potential, 
band edge of each TiC >2 layer and the energy level of Ti 
4s semi core orbital. The dotted lines are guide lines to 
show the potential flatness in the interior of TiC >2 slab. 
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FIG. S3. The step-like barrier at the pristine Au/TiC >2 inter¬ 
face. The solid curve shows the plane averaged electrostatic 
potential of the Au/Ti02 slab. The gray lines are the con¬ 
nection of the band edge of each Ti02 layer. The dots are 
the energy level of Ti 4s electrons in free standing and Au at¬ 
tached Ti02 slabs. The interface to the Au slab upshifts the 
energy level of Ti02 slab by 0.4 eV as a result of the interface 
electron polarization. 


NB-DOPANT INDUCED DIPOLE AND BAND 
BENDING 


The Nb-dopant induces a dipole field between the in¬ 
terface and the dopant layer. The magnitude was calcu¬ 
lated by integrating of Eq. 3 over the interface range and 
dopant layer. The induced dipole at the interface 
was found to be in reverse proportional to the dopant 
interface distance (do) as shown in Fig. [S4] 
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0 10 20 30 40 

d D (A) 

FIG. S4. The induced dipole at the interface {fAnd.) and 
around the Nb dopant (/Uj nl j) as a function of dopant inter¬ 
face distance (do). The dotted curve is fitted by = 

0.066 DA _1 /( d D - 1-08 A). 

The band bending is also presented by the shift of Ti 
4s semi core orbital energy. They show the same trend 
with the semiconductor band edges (Fig. IS5l) . 
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FIG. S5. The band bending in Ti02 varied with Nb dopant 
interface distance (do), (a) The layer localized conduction 
band minimum (CBM) referred to the Fermi energy, (b) The 
dopant induced energy shift of the Ti 4s semi core state. The 
vertical bars indicate the position of Nb-dopant in each calcu¬ 
lation. Both curves show that the TiCH bands shift up while 
approaching the interface. 

SCHOTTKY BARRIER HEIGHT PREDICTED 
BY MODELS 

In the analysis of SBH with image force (IF) correc¬ 
tion and deep level (DL) model, the dielectric constant 


Es was treated as a parameter. DL model with eg = 10 
shows good agreement with our first principles calcu¬ 
lations (Fig. 4). As TiC >2 has dielectric constant of 
4 000A = 97 and 58 in the c and a directions @], we also 
calculated the model predicted SBH reduction 
with dielectric constant of Es = 60 as shown in Fig. IS6l 

The uniform dopant Schottky model with IF correction 
(Eq. 4) is more sensitive to the dielectric constant as 
compared to localized DL model. With es = 60 the SBH 
reduction is always less than 0.2 eV in the considered 
density range. And the depletion layer width expands 
from 10 A (es = 10) to 25 A (eg = 60) under the dopant 
density of No = 1.3 x 10 21 cm -3 . 



d D (A) a D (10 14 cm 2 ) 


FIG. S6. The decrease of the SBH as a function of the dopant 
interface distance (a) and dopant areal concentration (b) us¬ 
ing model simulations with low and high dielectric constants 

(es)- 
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